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Abstract We have designed, fabricated and measured high-Q A/2 coplanar waveguide 
microwave resonators whose resonance frequency is made tunable with magnetic field 
by inserting a DC-SQUID array (including 1 or 7 SQUIDs) inside. Their tunability 
range is 30% of the zero field frequency. Their quality factor reaches up to 3x10^. We 
present a model based on thermal fiuctuations that accounts for the dependance of the 
quality factor with magnetic field. 
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1 Introduction 



On-chip high quality factor superconducting resonators have been extensively studied 
in the past years due to their potential interest for ultra-high sensitivity multi-pixel 
detection of radiation in the X-ray, optical and infrared domains [ll[2]. They consist of 
a strip line waveguide of well-defined length, coupled to measuring lines through input 
and output capacitors. The TEM modes they sustain have quality factors defined by 
the coupling capacitors and reaching in the best cases the 10^ range [2]. 

It has also been demonstrated recently that superconducting resonators provide 
very interesting tools for superconducting quantum bit circuits [4j. Indeed, a resonator 
can be used to measure the quantum state of a qubit [3][5ll6l[7] . Moreover, another res- 
onator may serve as a quantum bus and mediate a coherent interaction between qubits 
to which it is coupled. The use of resonators might thus lead to a scalable quantum 
computer architecture [5]. The coupling of two qubits mediated by a coplanar waveg- 
uide (CPW) resonator has already been demonstrated l8]|9]. In experiment each 
qubit needs to be tuned in and out of resonance with the resonator for the coupling to 
be effective. Reference [10] proposed an alternative solution that consists in tuning the 
resonator in and out of resonance with each qubit. Here we report on the measurement 
of high quality factor resonators whose frequency can be tuned. Measurements similar 



Quantronics Group, Service de Physique de I'Etat Condens (CNRS URA 2464), DRECAM, 
CEA-Saclay, 

91191 Gif-sur-Yvette, France 



2 





Fig. 1 a: Tunable resonator scheme : a DC SQUID array is inserted between two A/4 waveg- 
uides coupled to a 50 f2 measurement line through input and output capacitors Cc- b: Op- 
tical micrograph of a CPW niobium resonator, c: Typical coupling capacitor (design value : 
Cc = 27fF). d: Gap in the middle of the resonator, before SQUID patterning and deposition, 
e: Electron micrograph of an aluminum SQUID (sample A), fabricated using electron-beam 
lithography and double-angle evaporation, f: Electron micrograph of a 7-SQUID array (sample 
B). 



to ours have been reported by other groups on lumped element [TT] and distributed 
[T^ resonators. 



2 Tunable resonator with DC SQUID: model 

Our tunable resonators consist of A/2 coplanar waveguides with an array of A'^ DC- 
SQUIDs in series inserted in the middle of the central strip (see Fig. [1^). Each DC 
SQUID is a superconducting loop with self-inductance L; intersected by two nominally 
identical Josephson junctions of critical current Ico ; the loop is threaded by a magnetic 
flux The SQUID array behaves as a lumped non-linear inductance that depends on 
which allows to tune the resonance frequency. 

A A/2 CPW resonator without any SQUID consists of a transmission line of length 
/, capacitance and inductance per unit length C and £,, and characteristic impedance 
Zq = \/C/C. We consider here only the first resonance mode that happens when 
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Fig. 2 A DC SQUID with two junctions of critical current /^o E'nd loop inductance L; , biased 
by a magnetic flux ^ and by a current i, is equivalent to a lumped flux-dependent non-linear 
inductance Lj{4>, i) that can be decomposed in an inductance Ljo{<P) and a non-linear element 
SNL(^) in series. 



I — A/2 at a frequency cjr = tv/VLC, where L = CI and C = CI are the total 
inductance and capacitance of the resonator. The quality factor Q results from the 
coupling of the resonator to the i?o = 50 f2 measurement lines through the input and 
output capacitors Cc leading to 



from internal losses (Qint), and from possible inhomogeneous broadening mechanisms 
iQinh)- These combined mechanisms yield 

q-' = Qc' + q;^\ + q7„\- (2) 

As shown in Fig. (2] we model a SQUID as a non-linear inductance Lj{<l>,i) that 
depends on <1> and on the current i passing through it, so that the voltage across the 
SQUID is 



(3) 



All SQUID properties are periodic in <f with a period $q — h/2e, the supercon- 
ducting flux quantum. Introducing the reduced flux quantum ipQ = <1>q/2-k, the SQUID 
frustration / = Tr^/cfpi the effective critical current Icip) ~ 2Jcoj cos/| of the SQUID 
at zero loop inductance, and the parameter /3 = Lj/co/voi our calculation of Lj(<f, i) 
to first order in /3 and to second order in i/Idp) yields for / G ]— 7r/2, 7r/2[ 



Lj{-l>,i)^LjQ{$)+A{$)i' 



(4) 



with 



Equation U shows that the SQUID can be modelled as the series combination of a 
lumped inductance Ljq{(1>) and of a non-linear device SNL($) 13 (see Fig. [2)|. 

In the linear regime i <^ Ic{^) corresponding to low intra-cavity powers, one can 
neglect the non-linear term in Eq. [l] The N-SQUID array then simply behaves as a 
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lumped inductance NLjq($). The device works in that case as a tunable harmonic os- 
cillator. Introducing the ratio £(^) = Ljq{4>)/L between the total effective inductance 
of the SQUID and the resonator inductance, the frequency and quality factor are 



Qexti^) = Qcll + 'iNe{<P)]. (8) 

At larger peak current in the resonator i < Ic{^), the non-linear element SNL(^) 
has to be taken into account. The equation of motion of the oscillator acquires a cubic 
term, similar to that of a DufBng oscillator il4j. This leads to a small additional shift 
of the resonance frequency Su!o{E) proportional to the total electromagnetic energy E 
stored in the resonator. Retaining first order terms in e(#), we find 



Loom \nRo[l + 2Nsmj Slim ' ^' 

As shown by Eq. [T] a resonator including an array of A'^ SQUIDs of critical current 
NIco has approximately the same resonant frequency and same tunability range as 
a resonator including one SQUID of critical current Ico- However, an interesting ad- 
vantage of using an array is to obtain a linear regime that extends to larger currents, 
allowing measurements at larger powers and therefore higher signal-to-noise ratios. 



3 Sample fabrication 

The design and fabrication of our resonators closely followed Ref. [15] ■ The coupling 
capacitors were simulated using an electromagnetic solver. Test niobium resonators 
without any SQUIDs were first fabricated. They were patterned using optical lithog- 
raphy on a 200 nm thick niobium film sputtered on a high-resistivity (> 1000 Q cm) 
oxidized 2-inch silicon wafer. The niobium was etched away using either dry or wet 
etching. Dry etching was done in a plasma of pure SFg at a pressure of 0.3 mbar and 
at a power such that the self-bias voltage was 30 V and the etching rate 1.3nm/s. We 
observed that adding oxygen to the plasma gave consistently lower quality factors. 
Wet etching was done in a solution of HF, H2O, and FeCls having an etching rate 
of approximately Inm/s at room-temperature. A typical resonator and its coupling 
capacitor are shown in panels b and c of Fig. [T] Its 3.2 cm length yields a resonance 
frequency around 1.8 GHz. 

In addition to these test structures, some resonators had a gap in the middle (see 
Fig. [TJl) used in a later step to fabricate a SQUID array by e-beam lithography and 
double-angle aluminum deposition (see panels e and f in Fig. [1]). Before depositing 
the aluminum, the niobium surface was cleaned by argon ion-milling (dose < 10"'^^ 
neutralized 500 eV ions per square centimeter). The Nb/Al contact resistance was found 
to be in the ohm range, yielding tunnel junctions of negligible inductance compared to 
that of the SQUID. 
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Fig. 3 Experimental setup. The sample is thermally anchored at the mixing chamber of a 
dilution refrigerator with temperature 40 — 60mK. It is connected to a vector network analyzer 
(VNA) at room-temperature that measures the amplitude and phase of the S21 coefficient. 
The input line (top) is strongly attenuated (120 to 160 dB in total) with cold attenuators to 
protect the sample from external and thermal noise, and filtered above 2 GHz. The output line 
(bottom) includes a cryogenic amplifier with a 3K noise temperature and 3 cryogenic isolators. 
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Table 1 Summary of sample parameters. See text for definitions. 



4 Experimental setup 



The chips were glued on a TMMIO printed-circuit board (PCB). The input and output 
port of the resonator were wire-bonded to coplanar waveguides on the PCB, connected 
to coaxial cables via mini-SMP microwave launchers. The PCB was mounted in a 
copper box. The coefHcient (amplitude and phase) of the scattering matrix was 
measured as a function of frequency using a vector network analyzer. Test resonators 
were measured in a pumped ^He cryostat reaching temperatures of 1.3 K, with typical 
input power of —50 dBm and using room-temperature amplifiers. We measured internal 
quality factors up to 2x10^ with both etching methods. 

The tunable resonators were measured in a dilution refrigerator operated at 40 — 
60 mK, using the microwave setup shown in Fig. O The input line includes room- 
temperature and cold attenuators. The output line includes 3 cryogenic isolators, a 
cryogenic amplifier (from Berkshire) operated at 4K with a noise temperature of 3K, 
and additional room-temperature amplifiers. The attenuators and isolators protect the 
sample from external and thermal noise. This setup allowed to measure the sample 
with intra-cavity energies as small as a few photons in order to operate in the linear 
regime, corresponding to typical input powers of —140 dBm at the sample level. 




frequency (GHz) ^/^^ 

Fig. 4 (color online) a: Measured (thin line) amplitude (top) and phase (bottom) transmis- 
sion of sample A for ^ = and fit (bold line) yielding a quality factor Q = 3300. b: Measured 
resonance frequency of sample A (squares) as a function of applied magnetic flux and corre- 
sponding fit (full line) according to Eq. [7] 



5 Experimental results 

Two tunable resonators were measured: sample A has only one SQUID (see Fig. [1^) 
and large coupling capacitors (27 fF) so that its total quality factor is determined by 
Qc = 3.4 xlO^. Sample B has an array of 7 SQUIDs in series (see Fig. [Tf ) and smaller 
coupling capacitors (2fF) so that its quality factor is likely to be dominated by internal 
losses or inhomogeneous broadening. Relevant sample parameters are listed in table [1] 

A typical resonance curve, obtained with sample A at # = for an input power 
of —143 dBm corresponding to a mean photon number in the cavity n ~ 1.2, is shown 
in Fig. |4l The jS2i| curve was normalized to the maximum measured value. By fitting 
both the amplitude and the phase response of the resonator, we extract the resonance 
frequency and the quality factor Q. When the fiux through the SQUID is varied, the 
resonance frequency shifts periodically as shown in Fig. |4)d, as expected. 

The resonance frequency fo{$) and quality factor Q(^) are shown for both samples 
in Fig. [S] over one flux period. The /o(^) curves in panels a and c are fitted with 
Eq. [7] The agreement is good over the whole frequency range, which extends from 
1.3 to 1.75 GHz, yielding a tunability range of 30%. The small discrepancy observed 
for sample B might be due to a dispersion in the various SQUID loop areas that is 
not taken into account in our model. The parameters obtained by this procedure for 
both samples are shown in table [T] they are in good agreement with design values and 
test-structure measurements. 

The Q{$) dependance for both samples is shown in panels b and d of Fig. (5] Both 
samples show a similar behaviour: the quality factor depends weakly on <P when the 
flux is close to an integer number of fiux quanta, whereas it shows a pronounced dip 
around $o/2. 

The largest quality factors are 3.5x10'^ for sample A and 3x10^ for sample B. 
This difference is due to the different coupling capacitors. For sample A, the maximum 
quality factor is the same as measured on test resonators with similar capacitors and 
corresponds to the expected Qc for Cc ~ 27 fF. Therefore sample A quality factor is 
limited by the coupling to the 50 lines around integer values of (Pq. The situation 
is different for sample B: the measured value is one order of magnitude lower than 
both the quality factor Qc =6x10^ expected for Cc = 2 fF and the measured Q of 
test resonators with the same capacitors (see table [1} ■ This unexplained broadening 
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Fig. 5 (color online) a and c: Measured resonance frequency /o as a function o{'P/<Po (squares) 
for samples A and B, respectively, and fit according to Eq. [7] (solid line), b and d: Measured 
quality factor Q (disks) as a function of <P/<Po. The solid line is calculated according to the 
model (see text) for a temperature T = 60 mK. 



of the resonance in presence of a SQUID array might be due either to the presence 
of low-frequency noise in the sample, or to a dissipation source specifically associated 
with the SQUIDs. We note that flux-noise is not plausible since our measurements 
show no clear correlation with the sensitivity of the resonator to flux-noise. However, 
critical-current noise could produce such effect. Another possibility could be dielectric 
losses in the tunnel barriers. 

We now turn to the discussion of the dip in Q{$) observed around ^o/2. We 
attribute it to thermal noise. Indeed, as discussed in section[2l the resonance frequency 
depends on the energy stored in the resonator. At thermal equilibrium, fluctuations in 
the photon number translate into a fluctuation of the resonance frequency and cause an 
inhomogeneous broadening. At temperature T, the resonator stores an average energy 
given by Planck's formula i5 — hu}Q($)n,n = l/{exp[fta;o(^)/fer] — 1} being the average 
photon number. The photon number and energy fluctuations are — ri^ = n{n + 1) 
and 



' E +hojo('P)E. 



(10) 



The characteristic time of these energy fluctuations being given by the cavity damp- 
ing time Q/ujQ with Q 3> 1, a quasi-static analysis is valid and leads to an inhomoge- 
neous broadening Suji^h = |da;o/di?| \f5E^. Using Eq. [O] we get 



N 



x 



fo 



T:Ro[l + 2Nei'P)]j 81^ {<P) 



(11) 



The resulting quality factor is Q ^ = Qinh~^Qext' which is plotted as full curves in 
panels b and d of Fig. [5] for T = 60 mK. The agreement is good, although Eq. llll results 
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from a first-order expansion that is no longer valid in the close vicinity of /2- We have 
also observed that Q values significantly degrade around when the samples are 

heated, while remaining unchanged around integer numbers of ^g- These observations 
suggest that thermal noise is the dominant contribution to the drop of Q. Note that our 
model does not take into account flux-noise, which evidently contributes to Qinh ^nd 
could account for the residual discrepancy between experimental data and theoretical 
curves in panels b and d of Fig. [5] 

6 Conclusion 

We have designed and measured SQUID-based stripline resonators that can be tuned 
between 1.3 GHz and 1.75 GHz, with a maximum Q =3x10^ limited by an unknown 
mechanism. The quality factor degrades due to thermal noise around ^o/2. This lim- 
itation would be actually lifted with higher frequency resonators matching typical 
Josephson qubit frequencies. Their tunability range at high Q would then be wide 
enough to couple a large number of qubits. 
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